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ABSTRACT 


The soil profile is divided for the purpose of the study into layers 1 cm 
thick. The monovalent and divalent cations are treated separately in each 
layer. The flow of cations between the different layers is controlled by the 
movement of water and the ions in it. The values of the parameters are 
estimated using theoretical argumentation. The soil profile is assumed to be 
in equilibrium before the acid load is applied. The time constants of the 
cation exchange process are approximated. The equilibrium requirement and 
the time constants provide the information necessary for the estimation of 
the values of the parameters. The model is tested using lysimeter 
experiments with artificial acid rain. 


1. INTRODUCTION 


The effect of acid rain on soil takes place through weathering of minerals, 
microbial activity and cation exchange. The possibility of permanent change in 
the state of forest soils is most apparent in connection with the leaching of 
nutrients due to acidification of soil, where cation exchange is believed to 
play a major role. The leaching process can remove nutrients originally 
available for plants that are not recycled naturally, and another irreversible 
change can take place by the activation of certain poisonous processes such as 
release of aluminium compounds under acid conditions. 


In boreal coniferous forests the humus is relatively acid and the bulk of the 
decomposers are fungi that are rather insensitive to acidity (Trudgill 1977). 
This indicates that the impact of acid rain on micro-organisms is less 
significant than in more basic soil types. Secondly, the organic layer is thin 
and the buffer capacity of the soil is poor due to relatively low Ca content 
in the minerals. These factors indicate that in boreal coniferous forests the 
leaching of nutrients from the soil can be an important consequence of acid 
deposition. 
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Figure 1. Restriction of the system. 


The present paper introduces a simulation model of the acidification of soil, 
based on cation exchange and water percolation through the soil. The model 
produces an estimate of the change of soil acidity in time when the rate and 
acidity of precipitation are given. The work is part of a project dealing with 
the effects of air pollutants due to energy production on trees, carried out 
at Helsinki University. 


2. THE MODEL 
Fig. 1 illustrates the most important processes associated with soil nutrient 


status and acidity (cf. Trudgill 1977). Nutrients occur in the soil in four 
different states: in water as soluble ion (N3) adsorbed on the surface of 
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soil particles (Ni ), as chemical components of minerals (Ni ) and soil 
organic matter (N..). Hydrogen occurs either in soil water as Free soluble 
ions (H ) or adsorbed on the surface of soil particles (H.). The arrows in 
Fig. 1 describe nutrient and hydrogen flows between these° compartments and 
between the soil system and its environment. These flows only occur in the 
presence of water (W), which can either stay in the system or percolate 
through it from the surface to the lower soil layers. 


The aim of the study was to produce a model that can be utilized in analysing 
the effects of acid deposition on this system. The most important dynamic 
properties of the system are associated with the ion exchange process, 
provided that this is connected with the flow of water through the system. 
Therefore, the present preliminary version of the model only incorporates 
these two processes, ignoring decomposition of soil organic matter, weathering 
of minerals and nutrient uptake by plants. This simplification corresponds to 
the assumption that these three material flows approximately cancel each 
other. The state of the soil, X, is therefore defined as follows: 


m T 
X= (Hy Nps Nos Hye Wys Nogs W) 


H amount of hydrogen ions in soil solution 
amount of monovalent cations in soil solution 
N3 amount of divalent cations in soil solution 

amount of hydrogen ions adsorbed on soil particles 
amount of monovalent cations adsorbed on soil particles 
amount of divalent cations adsorbed on soil particles 
amount of water 


The amounts of ions are expressed in millimoles (mmol) and that of water in 
liters (1), both per unit volume of soil. 


2.1 lon exchange 


Ion exchange is treated analogously to the process occuring in a chemical 
buffer solution where ions with the same valence behave identically. 


Two types of cations are considered: N} and N represent monovalent 
cations sush as K and Na , while N dnd N Vapresent divalent ones 
such as Ca” and Mg’. The following réaction &fuations are obtained: (cf. 
Brustad et. al. 1979). 


+H, == tH (1) 
Nyt 2H, ee Ny, + 2H (2) 
2N) + Nog J 2N; +N (3) 


The error due to combining cations with the same valence is minor provided 
either that they behave similarly in the exchange process or that one of the 
ions clearly dominates the others. In leaching experiments by the Norwegian 
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SNSF project, the concentrations of cat and mgt changed similarly when 
water was percolating through a soil profile. The changes of K and Nw 
concentrations have ususally been similar to each other (Brustad et. al. 
1979). The experiments carried out in our project support these results. The 
positions of the ions with respect to nutrient uptake is different however, 
and if more detailed results are aimed at, different cations will have to be 
considered separately. 


The dynamic equations of the model are formed on the basis of Eqs. (1)-(3). We 
assume that the reaction rate is proportional to the amounts of the reacting 
ions. The differential equations formed on the basis of Eq. (1) for H and Ny 
are hence: 


dH 


dt S 9 HN, + ag HAN (4) 
aN, 
Mis ay HN, - ag HN (5) 


where a, and ap are constant parameters. 


Eqs. (2) and (3) are treated in a similar way, and the resulting derivatives 
are summed up to give a model for the complete system. The equations for H., 
and N. are obtained by changing the sign of the respective solubfe 


OA deri vatives . 


The following set of equations is obtained: 


ao, 
ge W~ ( aja + aglNs,)H +( aN) + cio Ha 
aN, 
et = -è oHa + ogNiNog)Hy + (aH + agNaN a )N a 


dN, = - 2 2 2 2 
ate C ogha + agNy,)Np + ( agh” + ash )Nog 


dH, __ dH (6) 
at? ~~ at 
dN, _ _ dN. 
atl? * - Gel 
dN, _ _ dN, 
ane? = - ate 


all state variables and parameters non-negative. 


296 


2.2 Transport of ions by water 


Consider the flow rate of water per unit surface area. In the top layers of 
soil horizontal movements of water are negligible, and from the point of view 
of leaching it is only the vertical dimension that matters. 


The model assumes that the water content of the soil never goes below field 
capacity, which corresponds to the watering experiment in the project and the 
situation in the field during a rainy period. It is further assumed that the 
surplus of water is tranported down the profile with a specific (constant) 
flow rate. The profile is divided into n horizontal layers, assuming a 
complete mixing of water and its solutants within each layer. The ions are 
hence transported down the soil profile whenever there is a downward flow of 
water. A summary of the water flow and ion transport model is presented in 
Figs. 2: 


letut yl 1 
y =s (W -W fo)e 


2 ep ye ye 2 
y^ =s (WE -W foe 


daft ast i 
y =s(W -W fee! 


n Myn n 
y =(W -W gc )e 


Figure 2. Description of the model of water percolation. Symbols: 
u,input (water and ions) 
ci vector of ion concentration in layer i 
y; output from layer i (water and ions) 
Wi amount of water in layer i 
W field capacity corresponding to layer i 
s bframeter 


2.3. Mathematical summary of the model 


Denote the time derivatives of H, Ny and No in Eq. (6) by DH, ON, and 
DN,, respectively. Associate a label~k (1,...,n) to each soil layer, the 
f16w of ions i with pater into the layer k by u` and the outward ion flow 
from the layer by y`. Combining the submodels for ion exchange and water 
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flow gives the following differential equations for hydrogen and exchangeable 
cations in layer k (the dot refers to time derivatives): 


ik = DH + uN - yf 
Eo ake uhh 
Ny = DNS + uNe = Yie 
fk = ak 

it, -ont 

ia < -0n3 

uk = ri x Skew 


There are seven parameters in the model, which have been determined combining 
a theoretical equilibrium analysis with some reaction rate estimates found in 
the literature. The differential equations have been solved numerically using 
the Runge-Kutta integration method. 


3, RESULTS AND DISCUSSION 

Simulations of the model have been carried out with initial hydrogen and 
nutrient ion concentrations corresponding to the values of VT-type pine 
forests, which is a common poor forest type in Finland (Table 1). The 
parameter values applied are listed in Table 2. 


Table 1. Ion concentrations in initial state in a 1 cm soil layer. 


Soil solution 


- amount of hydrogen ions H 1)0.47 mmol /m? 

- amount of monovalent cations Ny 1.30, ™ 

- amount of divalent cations No 0.48 " 
Adsorbed ions 

- hydrogen H 1200.00 " 

- monovalent cations Nia 50.00 " 

- divalent cations Nos 200.00 " 
Amount of water (field capacity) Wee 7:5 1/me 


1) Corresponds po pH 4.2 


298 


Table 2. Parameter values. 


a 0.729 1078 
al 0.110 1026 
as 0.595 1078 
a 0.380 1072 
ag 0.269 1078 
ap 0.758 10 

s 4.0 


ION CONCENTRATIONS IN LEACHATE 


H.Ny.N2 
0.74 (mmot /1) 
0.6 
0.5 
-H 
0.4 
—N: 
03 3 
0.2 Ng 
0.1 
0 i7) a pe 
0 10 20 30 


Number of showers 


ONE SHOWER 51 water/m2 


1200 s 
pH 3.0 


Figure 3. The simulted time course of ion concentrations in water flown 
through a 25 cm soil profile subject to,30 successive acid showers. 
The amount of water per shower is 5 1/m° over a 1200 s period and 
the pH value is 3.0. 


Fig. 3 shows the simulated time course of ion concentrations in water flown 
through a 25 cm soil profile Subject, to 30 successive acid showers. The amount 
of water (pH 3) per shower is 5 1/m* over a 1200 s period. In 20 showers the 
concentration of divalent cations reaches a maximum and starts to decline, 
while the concentrations of hydrogen and monovalent cations keep increasing 
till the end of the simulation. However, N, will achieve a similar peak if 
the simulation time is extended. This refldcts the replacement of nutrient 
cations by hydrogen ions on the surface of soil particles in cation exchange, 
which takes place until all the exchangeable nutrients have leached out. 
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ION CONCENTRATIONS IN SOIL SOLUTION 


Amount of ions Amount of ions 


(mmot/m2- cm) (mmol /m2- cm) 
H 
10 1.0 
Ng 
N 
pH 5.5 1 
05 0.5 pH 3.5 
HAN 
N2 
7 T m 
O 2000 4000 6000 8000 O 2000 4000 6000 8000 
Time (5) Time (s) 
Amont of ions H 
{mmol /m?- cm) 
254 
204 
154 
pH 2.5 
104 
54 
5 20000 40000 
Time (s) 


CONTINUOS RAIN, 4.2 ml/s 


Figure 4. The simulated time course of ion concentrations in the soil solution 
when the soil is exposed to continous acid precipitation. The input 
rate of water is 4.17 ml/s, applied over a period of 8000 s. The 
acidity of the precipitation is shown in the figure. 
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Whether or not it is the divalent cations that get exhausted first depends 
upon the subject soil, represented in the model by parameter values. 


Fig. 4 compares the development of ion concentrations in soils exposed to 
different degrees of continous acid precipitation. The input rate of water is 
4.17 ml/s, applied over a period of 8000 s. 


There is rather a distinct limiting acidity (pH 2.5 - 3.0) at which the buffer 
capacity of the system seems to collapse, resulting in a rapid increase in the 
amount of H ions in the solution and a nearly complete vanishing of 
nutrient cations. This is consistent with the experiments carried out in our 
project (Katainen and Heikkilä 1984). 


The relationship between the total amount of input hydrogen ions and the rate 
of input was analysed by comparing two simulations where identical amounts of 
hydrogen ions were applied at different rates. The result was that leaching 
was more efficient in the high rate simulation, indicating that the cation 
exchange process behaves exponentially with respect to H concentration. 

In summary, the model seems to function satisfactorily and hence to provide an 
applicable tool for further analyses. The most significant simplification was 
to ignore nutrient uptake by plants, and the results should be analysed 
bearing this in mind. The model apparently overestimates cation leaching 
during periods of high root activity, when the surplus of nutrients released 
by acidification is taken up by the roots. Considering the soil to be in field 
capacity has a similar effect - during dryer periods there is little 
throughflow of water and hence little output of cations. On the other hand, 
when the soil is dry the concentrations of elements are higher and the 
exchange process is more intensive, increasing the rate of change due to 
acidification. If there is a potential for the weathering of aluminium 
compounds, for instance, the model predicts it will be higher during dry 
periods. 


It is fairly easy to include a simple treatment of nutrient uptake and 
changing soil water content in the model, which will give more information of 
the significance of acid precipitation for forest soils. In order to 
incorporate tree growth, however, one needs a far more complicated approach, 
The point where the lack of knowledge seems to condence is the dynamics of 
root functioning in relation to soil characteristics. However, since this 
approach is restricted to mineral nutrients only, it is best utilizable trying 
to predict the rate of soit acidification as such. This is interesting for 
instance from the point of view of the poisonous reactions caused by the 
release of aluminium which only occurs under highly acid conditions. If the 
effects on forest flora are to be thoroughly examined, nitrogen and hence soil 
organic processes should be incorporated. 
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